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Abstract 

A knowledge of the acoustic energy emission 
of each blade row of a turbomachine is useful for 
estimating the overall noise level of the machine 
and for determining its discrete frequency noise 
content. Because of the close spacing between the 
rotor and stator of a compressor stage, the strong 
aerodynamic interactions between them have to be 
included in obtaining the resultant flow field. 

This paper outlines a three-dimensional theory for 
determining the discrete-frequency noise content 
of an axial compressor consisting of a rotor and a 
stator each with a finite number of blades. The 
lifting surface theory and the linearized equation 
of an ideal, nonsteady compressible fluid motion 
are used for thin blades of arbitrary cross sec- 
tion. The combined pressure field at a point of 
the fluid is constructed by linear addition of the 
rotor and stator solutions together with an inter- 
ference factor obtained by matching them for net 
zero vorticity behind the stage. The rotor solution 
is obtained as a Fourier-Bessel series; the stator 
solution is expressed as a Fourier-Laguerre series. 
The coefficients of the series and the interaction 
factor are determined as the eigenvector of a set 
of algebraic equations in matrix form whose ele- 
ments comprise the sum-integrals of the rotor and 
stator eigenfunctions. The resultant pressure 
field of the stage is the sum of the individual 
perturbation pressures in the presence of the 
interaction effects, expressions for which are 
given N herein. 

1, Introduction 

When the rotor is rotating in a uniform and 
steady incoming free stream, a point behind the 
rotor experiences periodicity in the pressure and 
velocity characteristics, which are propagated as 
acoustic perturbations. The relative motion 
between the stator blades and the rotor efflux 
introduces periodicity in the inflow to the stator. 
Since the rotor efflux consists of both the heli- 
cal vortex sheets and the viscous wakes of the 
rotor blades, the aerodynamic forces and moments 
on the stator blades and the outlet velocities and 
pressure of the stator are affected by the flow 
periodicity. Under subsonic flow conditions, peri- 
odicity in the stator flow field also affects the 
rotor. Reflection, transmission of the acoustic 
perturbation from the blade surface, and interac- 
tion with the duct resonance characteristics also 
have a strong effect on the resulting acoustic 
field at a point. 

Kemp and Sears** 2 studied the unsteady flow 
field of two rows of a two-dimensional thin airfoil 
cascade in incompressible flow although they did not 
dwell on the aeroacoustic aspects. The aeroacoustic 
field of two-dimensional single cascades was studied 
by several workers (Lane and Friedman 3 , Carta 4 , 

Carta and Fanti 5 , Hetherington 6 , Parker 7-9 , Tyler 
and Sofrin 10 * 11 , Smith 12 , Whitehead 13 , Morfey 14-16 , 
Mani 17 * 13 and Mani and Horvay 19 ). Kaji and 
Okazaki 20-22 studied the aerodynamic interaction 


between the two rows of infinite airfoil cascades 
and observed the discrete-frequency circumferential 
modes and their duct resonance characteristics. 
Slutsky 23 solved the discrete-frequency noise gen- 
eration due to rotor forcing excitation and rotor- 
stator interaction by using a Kemp-Sears mechanism 
with bound vortices on the rotor and the associ- 
ated three-dimensional velocity potential; he ob- 
tained the spin wave characteristics, Namba 24 * 25 
proposed the lifting surface theory for a single 
blade row in a compressible three-dimensional flow 
to calculate the unsteady blade forces and their . 
acoustic power generation. He found that the sound 
power generated was decreased by increasing the ra- 
dial nonuniformity of a sinusoidal external distur- 
bance in the circumferential direction. Mariano 23 
considered the effect of reflection of the sound 
field from a plane boundary for a uniform distribu- 
tion of sound sources. 

In an earlier paper 27 , the author outlined a 
lifting surface theory for a turbomachine stage in 
a three-dimensional incompressible flow that uses a 
distribution of sources and vortices to represent 
the blades of the two rows: the rotor and the 

stator. Another paper 23 treats the acceleration 
potential approach to the rotor-stator combination 
for compressible flow for both the aerodynamic 
loading and the acoustic characteristics. The con- 
cept of a rotor-stator interaction factor is intro- 
duced to isolate their mutual interaction effects, 
and axially attenuated or wavelike solutions are 
constructed for the general solution by using the 
Birnbaum series expansion as modulated by the 
eigenfunctions of the rotor and stator. Thus the 
composite solution of the rotor and stator is ex- 
pressed through Fourier-Bessel and Fourier-Laguerre 
series expansions that simultaneously satisfy the 
surface boundary conditions on both the blade rows 
and arbitrarily match to give zero net resultant 
vorticity in the wake of the stage. The procedure 
followed in this paper for two blade rows parallels 
closely that of Namba 25 for one blade row. An ad- 
ditional condition is required, however, to quan- 
tify the magnitude of the interference between the 
two blade rows by matching them optimally. In 
Section 2 the Green's function is written for the 
linearized differential equations governing the 
perturbation pressure field of a pulsating unit 
pressure monopole placed at a point on a blade in a 
double Fourier-Bessel and Fourier-Laguerre series 
for the rotor and stator, respectively. From the 
monopole solution the pulsating unit pressure 
dipole solution is obtained in Section 3 as the 
derivative normal to the blade corresponding to the 
surface pressure distribution. The chordwise pres- 
sure distribution on the blade is represented in 
Section 4 as a Birnbaum series modulated by the 
pressure eigenfunctions of the rotor and stator 
corresponding to the camber and thickness effects. 
The combined rotor-stator field is matched to give 
a zero net vorticity condition behind the stage 
(Section 5), and the coefficients of the Birnbaum 
expansion are determined by satisfying the flow 
tangency condition on the blade surface (Section 6) 
followed by a short discussion in Section 7 of the 
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combined pressure field. Because of space limita- 
tions, only the highlights of this work are pre- 
sented herein. No attempt has been made here to 
consider the reflection and absorption of waves or 
the natural duct resonance characteristics. Numer- 
ical results of the present theory will be reported 
separately. 

The author is deeply grateful for the keen 
interest evinced by Dr. Carl F. Lorenzo during the 
formulation of the present research program and for 
his constant encouragement throughout its execu- 
tion. Thanks are also due to Dr. Anatole Kurkov 
for his manifold discussions and useful suggestions 
on the various parts of the research and to Dr. 

James A. Pennline for his useful suggestions on 
some aspects of the solution of the Laguerre dif- 
ferential equation. The Aerodynamics and Engine 
Systems Division of the NASA Lewis Research Center 
provided the necessary atmosphere for the present 
work for which the National Research Council pro- 
vided the financial support. 

2 . Eigensolutions of Pulsating Pressure Pole 
on Rotor and Stator 

We shall consider a stage of an axial com- 
pressor (Fig. 1) with the rotor followed by a 
stator located centrally in an infinite annular 
duct with an outer diameter such that the clearance 
between the blade tip and the casing is uniform and 
negligibly small. We shall consider a compressible 
nonviscous fluid and assume a finite number of 
rigid blades Z r and Z s in the rotor and stator, 
respecti vely. Assuming that the compressor stage is 
lightly loaded so that the perturbation pressures 
are small compared with the free-stream static pres- 
sure poo, the linearized differential equations for 
the perturbation pressures p r and p s of the 
rotor and stator are given by 


velocity diagram (Fig. 2) so that we can write the 
components of inflow to the rotor and stator as 


( 2 ) 


(0, V r , W r ) = (0, nr, W fl ) 

(0, V s , W s ) = (0, nr - W a tan a 2r , W fl ) 

where the axial velocity W a is assumed to be 
constant throughout the stage. 


Equations (1) can be nondimensional ized, with 
all lengths and time expressed in terms of the 
rotor tip diameter and characteristic time, and 
written as 
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in cylindrical polar coordinates where (0, V r , W r ) 
and (0, V s , W s ) are assumed to be free-stream 
velocity components along the radial, peripheral, 
and axial directions for the rotor and stator, 
respectively, and a* is the free-stream speed of 
sound. Because the stator is situated downstream of 
the rotor in a region of rapid change, its inflow 
components (0, V s , W s ) cannot be defined a priori. 
However, we assume that the stator inlet conditions 
correspond approximately to that obtained from the 
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Equations (3) express the pressure perturbation 
field of a pulsating unit pressure pole situated on 
a rotor or stator blade at the vector point or 
with a pulsation frequency w- or (o $ , respec- 
tively. We may write the general solutions of 
Eqs . (3) as in Namba 25 

P r ( r i> p r ) = ~2 ^ f P r* r l ,p r ; W 
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x exp i[k r (e - 9p ) + a p (z x - c r )]da p (5a) 
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where we have defined 
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where a r and a s are the rotor and stator axial 
wave numbers; k r and k s (k r ,k $ = 0,+l,+2, . . . ,+«) 
are the pair of integer separation constants; and 
P r and P s are the unknown pressure functions. 

We shall express the delta functions in 
Eqs. (3) in a combined Fourier series - Fourier 
integral form as 
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Combining Eqs. (3), (5), (6), and (7) and equating 
the corresponding terms on the two sides, we obtain 
the nonhomogeneous second-order ordinary differen- 
tial equations for the radial pressure functions 
P r and P s , given by 
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An examination of Eqs. (3) or (8) shows that when 
x = 0 as in the case of an incompressible fluid or 
when the exit rotor blade angle a£ r =0 in a com- 
pressible fluid, both the equations reduce to the 
Bessel type. But, for all other cases, x cannot 
be neglected since k s extends over an infinite 
range of values, and it is essential to consider x 
in a formal solution of Eq. (8b). Thus Eq. (8a) is 
the nonhomogeneous Bessel differential equation and 
Eq. (8b) is the classical Whittaker nonhomogeneous 
differential equation, for which the complementary 
solutions can be written in the form 
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in which I ( x ) and K(x) represent the modified 
Bessel and Hankel functions, respectively, and 
W (x) is the Whittaker function. To write 
tne solution of the complete nonhomogeneous equa- 
tions, it is convenient to express the solution of 
the rotor function P^, in an infinite ^gp^nsion 

in terms of the normalized functions a r for 

(k ) 1 

for Eq. (10a). The function $7 r' is defined in 

the appendix. The solution for the stator function 
P s is more conveniently written in terms of the 
set of generalized Laguerre functions that satisfy 
certain orthogonality relations than directly in 
terms of the Whittaker function, to which they are 
related by the expressions 
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Thus we shall obtain the solutions of the nonhomo- 
geneous Eqs. (8) for p r and p $ in the form 
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where g^r' and hj^s^ are coefficients of the 
expansion to be determined. The eiqenf unctions 
r 1 and s ; are combinations of Bessel and 
Laguerre functions, respectively, that satisfy 
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in which the independent variables r £ S and q^s 
are defined by 
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For the Laguerre equation (Eq. 10a), eigensol utions 
exist if 31 » 32 
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is an integer, and we shall put 
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satisfy the orthogonality relations within tbe 
domain of rj, for each blade row such that 28 
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The Dirac delta functions a(ri - p r ) and 
6 ( r i - p s ) can be expanded in terms of the corre- 
sponding eigenfunctions such that 



The pressure functions P r and P s must satisfy 
the end conditions of zero radial flow at both the 
hub and the tip for which the radial pressure grad- 
ients must vanish at both these points and we have 
the relations 
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The perturbation pressure amplitudes p r and 
p s in Eqs. (5) involve summation over k r and k s 
and a contour integration over a r and a s with 
the contour so chosen as to yield bounded or expo- 
nentially decaying terms in the axial direction. 
Because of the cumbersome and lengthy nature of the 
calculations, only the final expressions for sub- 
sonic axial flow are given herein as 
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3. Representation of the Lifting Surface 

The lifting surface of the rotor and stator 
blades can be regarded as composed of a distribu- 
tion of pulsating pressure dipoles with constant 
interpole phase angle factors o r and a s and 
with the dipole axes normal to the local blade sur- 
face. The position of the unit dipole on the blade 
has the coordinates [p r , ^ + <p r + m r (a r + 6*0; 

z r + £r^i. [ p s> 9 S + + Ps^s + z s + ? s-l on 

the (m r ) tn rotor and (m s ) tn stator blade, respec- 
tively, in which <p r and <p s are the correspond- 
ing offset angles of the first blade (Fig. 3). The 
dipole field can be obtained from the pole field by 
differentiating along the local unit normal to the 
surface. The field of the Z r and Z s unit pres- 
sure dipoles on the blades of the rotor and stator 
can be obtained by summing over m r and m s . If 
Ap r and Ap s be the net upward-acting pressure at 
the points p^( P r ,¥ r ,c r ) and p" s (p s , 9 s ,c s ) , we can 
obtain the resulting pressure field of the rotor 
and stator each in isolation by integrating the 
dipole field of the intensity Ap r or Ap s and 
integrating over the entire surface of each blade 
with respect to p and For subsonic axial 
flow, we have the rotor and stator pulsating dipole 
pressure field given by 
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where the functions and y* are given in the 
appendix. 

4. Effect of Camber and Thickness 

The pressure loading Ap r and Ap s on 
the blades in isolation is caused by the combined 
effects of blade thickness, camber, and incidence 
and varies over the blade surface. It can be ex- 
pressed by combining it with the rotor and stator 
eigenfunctions V r ' and Vs' as 
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where (F^,Fp) and (Gj^Gp) are the chordwise varia- 
tions for the rotor and stator due to the thickness 
and camber, respectively, in terms of the chordwise 
coordinates w r and w s related to the locally 
orthogonal coordinates (yVr) and (ys» z s) 
defined by 
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The orthogonal coordinates (r W) and 
(r,ys,Zg) are related to the Cartesian coordi- 
nates (X, Y, Z) by rotation through the angles 
(a r ,a s ) and U r ,9s) described in Eq. (26) 
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The axial positions of the leading and trail- 
ing edges of the blades are given by 

z rl = z r0 - C R cos V z r2 = z r0 + C R cos “r 

(27a) 

z sl = 2 s0 - C R cos “s’ z s2 = z s0 + C S C0S “s 

(27b) 


where the coefficients Aj(, k r), B^r), C^s), and 
Djf/s) can be determined by making the blade sur- 
faces as stream surfaces. The series in Eqs. (29) 
satisfy the Kutta condition at the blade trailing 
edge and have the usual singularity at the leading 
edge. It can be observed from Eqs. (24) that 
thickness and camber effects are inseparable and 
are interwoven. Note that the thickness effect 
corresponds to the monopole solution and can be 
obtained directly. 


where (\|> r , 4 > $ ) and (c r ,c s ) are related to w r by 
the equations 



2 irm _i f sin a cos “ i 

+ cos \ j 


? r = Z r,0 + 2c R C0S “r cos “r 


? s = Z s,0 + 2c S cos “s cos “s 


(28) 


5. Boundary Conditions 

By defining the mean camberline and thickness 
shapes of the blades through the equations 
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The upper and lower surface blade profiles are 
given by 
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The offset angle of the first blade on the 
rotor can be written as = 1 * 0 * 1 * where 
the offset angle 9 / of the first blade on the 
stator can be set equal to zero. The terms 
(2irm r /Z r ) and (2irm r /Z s ) represent the polar angles 
of the (m r ) th rotor and (mg)* h stator blades, 
respectively, at the midpoint of the blade chord. 
The last terms in and give the increment 
in the polar angles for points (y‘,z') away from 
the midchord point at any radius. We shall repre- 
sent the functions (Fi,F 2 ) and (G^ ,G 2 ) in an infi- 
nite trigonometric series expansion corresponding 
to the Glauert/Birnbaum series as 
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If U r and U s be the resultants of the free- 
stream and induced velocity perturbations of the 
rotor-stator flow field, we can represent the sur- 
face flow tangency condition on the upper and lower 
blade surfaces by the equations 



(32b) 


In keeping with the thin airfoil theory postulates, 
the boundary conditions are satisfied on the blade 
chord. 

To obtain the induced velocity perturbations of 
the rotor and stator, it is convenient to introduce 
a local intrinsic helical orthogonal coordinate 
system (Fig. 4) for each based on the corresponding 
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undisturbed streamline with the helix angles e hi 
and e^2 defined by 

R* 

tan e hi = ^ 

(33) 


SI 


1 0 0 


v r 

Sy‘ 


0 sin a s COS a $ 


v 0 - tan “2 + ^ 

l Sz 1 


0 -COS a sin a _ 


V + 1 


s s 


z 


tan e h2 = 


R* 

(r 1 - M 2 /R) 


(37b) 


The helical intrinsic coordinates (p r , tj, s^,) 
and (p s , T 2 , S 2 ) of the rotor and stator are 
related to the cylindrical coordinates by the 
matrix transformations of differential lengths as 



The equations of motion in the intrinsic coordi- 
nates can be written as 


oj'' 

< 

-s 

CD 

< 

1 — 1 

CD 

< 

Q 

t— 

1 

/aP 

3 P 

\ 'N 

3 P \ 

[ r 

r 

r 1 

yas-^ 5 3s 1 ’ y 

_ rM r 

V r l* 

3t 1* 

3S l/ 


> (35) 


/ 3V r2 

3 \2 
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1 

/ 3P s 

3P s 3 M 

l 8S 2 ’ 

3S 2 5 

3S 2 / , 

rM s 

V r 2* 

3t 2’ 3S 2/ 


where M r and M s are the resultant free-stream 
Mach numbers of the rotor and stator, defined by 



J 


The velocity perturbations (v r i, v p i, v T i) and (v r 2, 
v p2> v t 2^ the ro ’ tor anc * the stator separately can 
be obtained by integrating Eqs. (35) along a helical 
streamline by putting e,. = (0 - z,/ R*) and inte- 
grating over zj_. We tdius^ obtain tne resultant 
free-stream velocity UJ.,lT s from the equations 


U R1 


1 0 

0 


v r 

U Ry' 

= 

0 sin 

COS a 

r 


v e + r l /R * 

U Rz' 


0 -COS a r 

sin a 

r 


v z + 1 


(37a) 


where the total perturbation velocity of the rotor- 
stator combination (v r , v Q , v 2 ) is expressed as 



The first two vectors on the right of Eq. (38) are 
the perturbation velocities of the rotor and 
stator, each in isolation, and the last vector is 
the interference velocity due to the simultaneous 
presence in the compressor stage of the two rows in 
close proximity. We shall assume that the inter- 
ference velocity vector is a linear combination of 
the perturbation velocities due to the rotor and 
stator, so that Eq. (38) can be written 



This is the simplest assumption that can be made 
for the interference effect and permits us to esti- 
mate the magnitude of the interaction. By defining 
the parameters 


®1 " 9 hl V e 3 = 0 hl " “s’ e l 
e 2 = e h2 - “r ; e 4 = e h2 " “s 


1/(1 + e);) 


(40) 


and combining Eqs. (32) and (40), we can obtain the 
perturbations v a j, v T i, v a 2 , and v T 2 in the form 



(41) 
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where the matrices [AA] and [BB] are given by 


AA = 


- (sin 

CD 

+ 

T U1 

cos 

ep - (sin 

e 2 

+ 

T U1 
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e 2 )(cos 

9 1 - T U1 
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e-^Mcos 

e 2 

■ T U1 
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e 2 ) 

- (sin 

6 1 

+ 
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e 2 
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6^ ) ( COS 

e 2 

" \1 
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- (sin 

e 3 

+ 

t U2 

cos 

83 ) - (sin 

0 4 

+ 

t U2 
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84 ) (cos 

0 3 " t U2 
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83 ) (cos 

0 4 

" t U2 

sin 
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- (sin 

e 3 

+ 

t L2 
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03 ) - (sin 

e 4 

+ 

t L2 
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0 3 " t L2 
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e 4 

' t L2 
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f [ ( T u 1 si n « r + cos o p ) (rj/R*) + (t^cos « r - sin a r )] 
C ( t, iSin a + cos a ) (r,/Rj + (t. ,cos o„ - sin a )] 


z ;=o+ 


\ 


z r = 0" 


[( t u2 sin « s + cos a $ ) ( r^/R* - tan « 2 ) + (t u2 cos a $ - sin a $ )] 


zj.=0+ 




C(t l2 sin a $ + COS a $ ) (r^R* - tan a 2 ) + (x L2 COS a $ - sin a $ )] 

" 0 J 

(43) 


( 42 ) 


^1 ,m 

^1 ,m 




f 2,m 

^2 ,m 

^(k) 

m 

= e i 


^3 ,m 

g 3,m 

^(k) 

m 


f 4,m 

g 4,m 



*4 


The velocity perturbations v a i, v a 2 , v T j, and v T 2 
can be expressed in terms of the coefficients of 
the Glauert series expansions and written as the 
system of equations in matrix form 


(44) 


K = 0, 

The submatrices f , g, and <€ in Eq. (44) are 
functions of the position vector and are not given 
here because of space limitations but may be ob- 
tained from Ramachandra^. The coefficients 

(k \ m 
and ' (m = 0,1, 2,...,®) are defined in terms 

of A^ k *, B? k t C« k ', and as 

^ . «' k > * B< k >;^ k > - «J k > -2B< k > ; 

m - % 4 ^0 ” L 0 u 0 * 


aJf< k > , c{ k > - 2D< k >;«< k > , C (k » - 0 <k > 

1 1 0 m m m 

m ~ 2,3,... 

6. Rotor-Stator Matching and Final Solution 


(45) 


With each of the rotor and stator blades is 
associated a bound vortex of strength r r ,r s pro- 
ducing the local lift L = pVr per unit span 
according to the Kutta-Joukowski theorem. The 
strengths r r and r s of the bound vortex can be 
obtained from a chordwise integration of Ap r and 
Ap s as a function of r^ and z ^ as 


V» W a 1 + 


z r 2 

!■/ 


L rl 


T.p VL 
s °° a 


1 + 


(ft - tan "a)' ' J 


*p r dz r 


L r 2 


(46a) 


* P s d2 s 




(46b) 

We assume that the rotor and stator are matched at 
any operating condition when the net vorticity 
behind the stator is zero. If r r h and r s h are 
bound vortex strength of the rotor and stator at 
the blade hub, we express the matching condition by 
the relation 


Z r r rh + Z s r sh = 0 


(47) 


By combining Eqs. (24), (29), (45), and (46), we 
can express the matching condition (Eq. (46)) as 


^ 0 ( f 5,rfi a m k) + 95,n$i k) ) = 0 


k = 0,±1,*2 *- 

(48) 

The reader is referred to Ramachandra^ for the 
expressions of f 5 and % in Eq. (48). A composite 
equation can be written in matrix form for the 
unknowns and by combining Eqs. (44) 

and (48) as 


k ) * gfi k ^ = o 


(49) 
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8. Discussions and Conclusions 


,(k) _ ^(kUk) . T 


where aM 

** = -“"m e l 

and is the matrix 


r (k) 


is the vector of unknowns 

m = 0,1,2,. .. 


(BO) 


f l,m 

g l ,m 

-«i 

f 2,m 

g 2,m 

-*2 

f 3,m 

g 3,m 

"*3 

^4,m 

g 4,m 


f 5,m 

g 5,m 

- 0 


(51) 


The eigenvectors of the matrix.^ 
give the coefficients 


i 


m = 0,1,2,. . . 

in Eq. (49) 
and S\ 

K ^ 0,±1,±2,.. .,*«) from which the pressure load- 
ing ap r and Ap s on the blades is obtained. 
Equations (45) are insufficient to determine the 
coefficients A, B, C, and D of the' Glauert 
series of Eq. (29) in order to separate the thick- 
ness and camber contributions. 

7. Combined Pressure Field 


It is now possible to express the pulsating 
pressure fields Pr and Pj of the rotor and 
stator from Eq. (23) inclusive of the interaction 
effects. The resultant pressure P at any point 
in the flow is the sum of the pressures Pr and 

PS 


P(r^,t^; k^,k s ,£; ay,u) s ) - Pp(r^,t^; k r ,t; ay) 

+ VW k s ,s,; "s } < 52 ) 

Assuming that the free-stream values of the stream- 
wise components of the perturbation velocities 
Vol,~ = v o2 ,«> = 0, we can write Eq. (52) in the 
form 

(4»B/Z r )P(r 1 .t; k r ,k $ ,i; uy,w s ) 


iw t o Vt, 

me r 1 - t^4$ J MZ s /Z r )e s 


j,£.,m=0 


where 


(53) 


(54) 


The terms (^ 1 ,^ 1 ), (%-^h and are > 

respectively, functions of r^, e, and z^ alone; 
#i, t[ 9 ^p, and are constants obtained in 
terms of trie geometry of the stage, the number of 
blades, and the flow conditions. Again, because 
of space limitations, the reader is referred to 
Ramachandra^B for detailed expressions of these 
quantities. 

The acoustic modes and their characteristics, 
including their attenuation and resonance proper- 
ties in ducts,/ can be obtained from Eq. (53). How- 
ever, these discussions will be given in separate 
papers elsewhere. 


A lifting surface theory has been formulated 
for the two rotor-stator blade rows of an axial 
compressor stage in terms of the acceleration po- 
tential. Since the acceleration potential equals 
the ratio of the perturbation pressure and the mean 
flow density, the method is equivalent to the deter- 
mination of the acoustic field of a rotor-stator 
combination in a duct of infinite length in terms of 
the flow Mach number, the geometry, and the oscil- 
lation frequency of a distribution of pulsating 
pressure dipoles through a lifting surface theory. 
The combined effects of blade thickness and camber 
have been included as well as that of incidence 
although it has not been possible to separate them. 
Curiously, it has been possible to separate their 
effects when the fluid is incompressible 2 '. The 
perturbation pressures designated in terms of their 
mode numbers (k r , k s , and a) form a superposition 
of an infinite number of acoustic modes consisting 
of the circumferential wave numbers kv* and k s and 
the radial mode numbers n. It has also been pos- 
sible to isolate the rotor-stator interaction 
effects by matching the blade rows to give zero 
vorticity downstream of the stator as in an ideal 
combination. This assumption may appear restric- 
tive in a multistage compressor. However, this 
can be relaxed and the same method can be used to 
determine the coefficient vector by using a finite 
prescribed wake vorticity condition in Eq. (47). 

From this general solution it is possible to obtain 
special cases of interest like that of a single or 
twin tandem actuator disk or single or twin two- 
dimensional cascade rows. This procedure could be 
extended to the study of counterrotating propel- 
lers, of course, by modifying the end conditions 
in Eq. (20) at the blade tip. This results in a 
different set of eigenvalues and the corresponding 
eigenfunctions. 


APPENDIX 


In this appendix expressions for the quan- 
tities mentioned in the paper are given without 
their detailed derivations. The functions I and K 
are, respectively, modified Bessel and Hankel func- 
tions of the first kind normalized over the range 

(k r > 

h < n < 1. The radial eigenfunction $ 
given in Eq. (Al) appears somewhat different from 
33 

the form used by McCune for the rotor. The 

(k ) (k ) 

radial eigenfunctions $ and y of the 
rotor and stator are given by 


(kJ 
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(k +1) (kJ 
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og 
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I (kr) (x (k ' r) r 1 ) + K ' (x r r x ) (Al) 
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where 


5 < k s> 


VI 2 / 
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(A3) 


x » -5-3/2 


(vj) 


The eigenvalues x(V) and x^ k s) are obtained as 
roots of the equations 


(k-1) (k+1) 

K r (x x ) + K (Xj) 
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where CF and SF are functions defined in 
Ramachandra^ and will not be given here because 
of their lengthy nature. 


(A5) 
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